Fano-Hopfield model and photonic band gaps for an arbitrary atomic lattice 



0^ 
O 
O 
(N 



Mauro Antczza and Yvan Castin 
Laboratoire Kastler Brossel, Ecole Normale Superieure, 
CNRS and UPMC, U rue Lhomond, 75231 Paris, France 
(Dated: July 22, 2009) 

We study the light dispersion relation in a periodic ensemble of atoms at fixed positions in the 
Fano-Hopfield model (the atomic dipole being modeled with harmonic oscillators). Compared to 
earlier works, we do not restrict to cubic lattices, and we do not regularize the theory by hand 
but we renormalize it in a systematic way using a Gaussian cut-olT in momentum space. Whereas 
no omnidirectional spectral gap is known for light in a Bravais atomic lattice, we find that, for a 
wide range of parameters, an omnidirectional gap occurs in a diamond atomic lattice, which may 
be realized in an experiment with ultra-cold atoms. The long-wavelength limit of the theory also 
provides a Lorentz-Lorenz (or Clausius-Mossotti) relation for an arbitrary lattice. 
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I. INTRODUCTION 

The determination of the spectrum of light in a pe- 
riodic ensemble of atoms is a fundamental problem that 
still raises intriguing questions. After the seminal work of 
Hopfield [H (see also Agranovich 0), based on the Fano 
oscillator model for the atomic dipole many theo- 
retical works have been performed Most of them, in- 
spired by the typical context of condensed matter physics 
considered in focus on the long- wavelength limit 
where Lorentz-Lorenz (or Clausius-Mossotti) type rela- 
tions may be derived d, @ . Recently, this problem was 
extended to the whole Brillouin zone in the case of cubic 
lattices [3, @] , which allows to address the presence or the 
absence of an omnidirectional spectral gap for light. 

This problem of light spectrum in atomic lattices is 
no longer a purely theoretical issue. Recent experiments 
with ultracold atoms, having led to the observation of a 
Mott phase with one atom per lattice site , have indeed 
opened up the possibility to investigate the propagation 
of light in an atomic lattice, taking advantage of the large 
variety of optical lattices that may be realized to trap the 
atoms , all this in the regime where the lattice spacing 
is of the order of the optical wavelength, so that a probing 
of the whole Brillouin zone can be envisaged. 

In the photonic crystals in solid state systems, made of 
dielectric spatially extended objects (rather than point- 
like atoms) , after the pioneering work of [Tl1 | for diamond 
lattices of dielectric spheres, many configurations are now 
known to lead to a spectral gap, with a variety of appli- 
cations to light trapping and guiding [l^, On the 
contrary, for atomic lattices, no omnidirectional spectral 
gap was found, neither in cubic atomic lattices 0] nor in 
several less symmetric Bravais lattices [l^ . 

Two factors determine the presence of an omnidirec- 
tional gap: crystal geometry and details of the light- 
matter scattering process, both can separately close a 
gap. Indeed, in photonic crystals materials, character- 
ized by a macroscopic modulation of the refractive in- 
dex, the same lattice geometry can lead or not to an 
omnidirectional gap depending on the modulation of the 



refractive index, as is the case for the simple cubic (sc) or 
body centered cubic (bee) lattices [l^. Differently from 
photonic crystals materials, in the physics of ultracold 
atoms, one can realize periodic structures with a single 
atom per site 0] rather than a macroscopic number, re- 
alizing an ideal crystalline structure with a variety of 
possible geometries [l^ . Atoms are scatterers character- 
ized by a strong resonant and point-like interaction with 
light, so that the features of the light propagation in an 
atomic lattice cannot be straightforwardly extrapolated 
from known results in solid state photonic crystals. Also 
in the atomic case, it is not possible, only by geometric 
considerations, to predict the presence of an omnidirec- 
tional photonic band gap, the details of the light-matter 
scattering process do matter Q. 

Here we develop, within the Fano-Hopfield model, a 
self-consistent theory for the elementary excitation spec- 
trum of the light-atom field in atomic periodic systems, 
which is valid not only for cubic symmetry atomic lattices 
3 but also for any Bravais lattice, and, even more, also 
for periodic non-Bravais lattices (i.e. for crystals with sev- 
eral atoms per primitive unit cell). Our theory includes 
the light polarizations degrees of freedom, and is based on 
the introduction of a Gaussian momentum cut-off which 
allows to eliminate all divergences (even in a periodic 
infinite system of atoms) by a systematic renormaliza- 
tion procedure. We then use our theory to address the 
existence of an omnidirectional spectral gap for light in 
atomic lattices. In particular, we show that the diamond 
atomic lattice, which is a non-Bravais lattice composed 
of two identical face-centered cubic (fee) atomic lattices 
deduced one from the other by a translation along the 
main diagonal of the cube, may support an omnidirec- 
tional gap for light. To our knowledge, this is the first 
example of the occurrence of such a gap in a periodic 
structure of point-like atomic scatterers. 

The paper is organized as follows. We present the 
model Hamiltonian in section [TTl We renormalize the 
model and obtain an implicit equation for the light spec- 
trum, first for Bravais lattices in section UTTl (with a com- 
parison to existing predictions for a fee lattice) , and then 
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for a general lattice in section IIVI We derive Lorentz- 
Lorenz relations in the long wavelength limit for an ar- 
bitrary lattice in section |Vl In section IVII we calculate 
the spectrum in a diamond atomic lattice, we discuss the 
existence of a gap and experimental issues such as the 
practical realization of a diamond atomic lattice. We 
conclude in section IVlTl 



II. THE MODEL 

We consider N atoms with a dipolar coupling to the 
electromagnetic field, each atom being modeled by har- 
monic oscillators in the spirit of the Fano-Hopfield model 
[U, 0] • The atoms have fixed positions on a Bravais lat- 
tice (in scction lTVl we will consider the even more general 
case of non-Bravais lattices). The i*^ atom is in lattice 
site Ri. The Hamiltonian may be written as 

i=l a k6l>e±k 

N 

- ^D,.]Ei(RO. (1) 

i=l 



and to the operator (r) 

E±(r)= J d^uE^{r-u)x{u;b), (3) 

which is the convolution of the transverse electric field 
operator 

E±(r) = E E [^1^^ e^ '' + h.c] (4) 

kGPeJ-k 

with a normalized cut-off function x(r; b), J (Pr x(r; b) = 
1, where £k = *[ftfcc/(2£oV)]^/^, and the length b is the 
cut-off parameter. Here x(r; b) is a real and even function 
of r. This convolution regularizes the theory by eliminat- 
ing ultraviolet divergences; it makes the model Hamilto- 
nian H well defined. From ^ and ^ one has that 

l^^(^) = E E [^fc^ + h.c] x(k; 6), (5) 

where the Fourier transform of x(r; b) with respect to r, 
X(k;6)= J d?r e-^''-^ xiv.b), (6) 
is a real and even function of k. 



The first two terms in ([T|) correspond to the uncoupled 
atomic and radiation contributions respectively, while the 
last term is the dipolar coupling operator. The operators 
bi^a and b\ ^ respectively destroy and create an atomic 
excitation for the atomic dipole i along spatial direc- 
tion a. They obey the bosonic commutation relations 
[bi,a,b]j^f3] = Sij6a/3 and = 0. In (P), ujb is the 

bare atomic resonance frequency, the sum J2a over the 
three directions of space x, y, z accounts for the three 
spatial components of the dipoles. 

The photon annihilation and creation operators also 
obey usual bosonic commutation relations such as 

[akejO-k/g'] = i^ee'^kk' and [ake,ak'E'] = 0, where e and 
k are the photon polarization and wavevector. c is the 
velocity of light in vacuum. We assume a quantization 
volume V which includes N = Ylj atoms and cor- 
responds to periodic boundary conditions of the field, 
with a period M^e^ along each direction 7 S {1,2,3} 
of the Bravais lattice, with G N and is one of 
the three basis vectors of the Bravais lattice in direct 
space. As the consequence, the allowed wavevectors k 
for the electromagnetic field belong to the discrete set 
I? = {k I k = ^^(m^/M-),) e^, Vm-y g Z}, where the 
e^'s are basis vectors of the reciprocal lattice, such that 
• e^i — lixb^^i. The atom-electromagnetic field cou- 
pling term involves the dipole operator of the i*'^ 
atom, whose component along direction a 



D, 



(2) 



III. LIGHT DISPERSION RELATION IN A 
BRAVAIS LATTICE 

Starting from the Hamiltonian ([1]) it is possible (as 
done in Appendix [X| to derive the equations of motion 
for the matter and light fields, and, after a linear trans- 
formation and the use of the Bloch theorem, furthermore 
leaving out the so-called free-field solutions (see below), 
one obtains that the atom-light elementary frequency 
spectrum uj, in terms of bare quantities, is given by the 
solutions of the equation 

deti\/(&) = (7) 

where the 3x3 symmetric matrix M is given by 



E 

KeRL 



{ujjcf - K 



K' K' 
" ''x^(K';&), (8) 



12 



q is the Bloch wavevector in the first Brillouin zone, K 
is a vector of the reciprocal lattice (RL), K' = K — q. 



2 _ ^dgLUB 



heoVL ' 



(9) 



is proportional to the bare atomic dipole moment d^. 



is the squared bare plasmon frequency, Vl is the volume 
of the primitive unit cell of the lattice L in real space. 
We have taken the infinite quantization volume limit, so 
that V ^ M.^, and the Bloch vector q may assume any 
value in the first Brillouin zone. 
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Now, we specify the cut-off function x(k; b) and renor- 
malize the theory by expressing the equations in terms of 
physical reahstic quantities ujo and d, instead of the bare 
ones ujb and d^. Inspired by p^ . we choose a Gaussian 
cut-off function [l^, [l^ 



-2(k;6)=e-'=- 



X 



b > 0, 



(10) 



which leads to a rapidly convergent sum in ([H]). If 6 ^ 
0, which is the limit we are interested in, the sum in 
([5]) diverges. The key idea is then that Eq.® has been 
expressed in terms of hare quantities. In what follows, we 
renormalizc the matrix Al^fs by collecting the divergent 
terms and the bare quantities in new physical quantities. 



Renormalization of the theory 

To this end, it is useful to introduce the quantity 
Aq/3(6) defined as: 



^apib) 



OafS ^ X 

DTT 

+ 2{ub/c)^ 



Erfi (Lub/c) 



27ri/2(w6/c)3 



= (wh/c)2 



KeRL 



(11) 



where K' = K — q and Erfi (x) — 2n^^^^ dy expf 
is the imaginary error function. It has been shown in [1^ 
that Aq,_3(6), which does not depend on bare quantities, 
has a finite limit Aq^(O), as 6 ^ [see Eq. p^ ]. It is 
then possible to write the matrix Maf}{b) in the limit of 
& ^ as: 



heo 



Vl 



Sal: 



5ap {Uj/cf 



127r3/263 



67r3/26 



(12) 



where we have used Erfi(cjfo/c) = 0{b). The key point is 
that, although the last two terms in the square bracket 
[. . .] of (|12p diverge as & 0, they diverge in a way 
not explicitly depending on the Bloch wavevector q [19| . 
In simple words, considering the fact that the first term 
in the right-hand side of Eq. (fT2|) is (a>^ — ijj'^)6afi'- The 
term inside the [. . .] diverging as 1/b^ may be summed to 
i^B^aiS, resulting in a renormalization of cij^; the term 
inside the [. . .] diverging as u'^/b may be summed to 
—uj'^Sap, requiring also a renormalization of the atomic 
dipole moment ds appearing in the overall factor to the 
left of [...]. The fact that two parameters are renormal- 
ized in the Hamiltonian H is reminiscent of the renor- 
malization procedure of QED, where both the electron 
charge and mass arc renormalized (20| . 



In practice, to apply the renormalization procedure, 
we pull out in the coefficient of the — term: 



Ma,3ib) = 1 + 



dgUJB 



37r3/2feQc2& 



+ Oib), (13) 



where, by construction, the 3x3 symmetric matrix M 
has a simple expression in terms of the renormalized res- 
onance frequency ujo and the renormalized dipole d: 



with the renormalized plasmon frequency 
2 2d^(j;o 



(14) 



(15) 



Identification of (fT3|l with (fT2|) leads to the explicit ex- 
pressions for the two renormalized quantities 



Wo 



67r3/2feoa;B63 



Sn^/^hsoc^b 



,2 d| 
d = Wo — 

UJB 



1 



67r3/2?ieQ(jj3f)3 



-1 
(16) 

(17) 



Finally, neglecting 0{b) in we obtain that the atom- 
light elementary excitations for a generic Bravais lattice 
are the solution of the equation 



det] 







(18) 



where M is given by (fT4|) and does not depend on the 
cut-off b. In practical calculations of the matrix elements 
of M, one uses the fact that AapiO) = limft^n A„/^(&) is 
related to A„^(6) of Eg. pTj) by the useful expression 



A„^(0)=A„^(6) e("^/^)' 



(19) 



valid for all values of 6 <C amin, where Omin is the min- 
imum distance between two atoms. In what follows, it 
will be convenient to use as a parameter, rather than the 
plasmon frequency, the free space spontaneous emission 
rate of a single atom. 



r = 



d2fc3 



(20) 



37reo^ 

where fco = ujq/c is the resonant wavevector. 

Free-field solutions 

An important note is that, as already mentioned, the 
solutions of ([7]) do not exhaust the spectrum. One should 
add the so-called free-field solutions located on the free 
photon dispersion relations and corresponding to a free 
electric field of frequency uj (a solution of Maxwell's equa- 
tions in the absence of matter) that vanishes exactly in 
all the atomic positions. This was discussed in Q and 
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also appears in the calculations of Appendix 1X1 In view 
of the discussion of spectral gaps to come, it is useful to 
keep in mind the following lower bound on the frequency 
tjfrcc of all possible free- field solutions [l^ : 



Wfr, 



> inf ^. 

KeRL* 2 



(21) 



This inequality holds for an arbitrary (non necessarily 
Bravais) lattice, and it is saturated for a Bravais lattice, 
as one can show using the two- mode ansatz of Q. 



Spectrum for a fee lattiee 

As a straightforward application of our approach, we 
calculate the light dispersion relation in the fee lattice. 
There are two existing predictions, 0] and [1] , that reach 
opposite conclusions concerning the existence of an om- 
nidirectional gap. We have calculated the dispersion re- 
lation for exactly the same parameters as these two ref- 
erences, see Fig in for and FiglH for Q respectively, 
which allows a direct comparison. Our theory disagrees 
with the result of ; it quantitatively agrees with the one 
of [1] provided one replaces in the final result of (but 
of course not in the Lagrangian of the bare values 
of the atomic and plasmon frequencies by their renor- 
malized values. We conclude, as in that there is no 
spectral gap for light in a fee atomic lattice. 



IV. LIGHT DISPERSION RELATION IN 
ARBITRARY ATOMIC CRYSTALS 

Let us consider the generic case of a crystal which is 
not a Bravais lattice, that is it has more than one atom 
per primitive unit cell; it can be seen as the composition 
of P translated copies of the same Bravais lattice. The 
crystal is then the periodic repetition of an elementary 
base of P atoms. In this case it is possible to generalize 
the equations obtained for a Bravais lattice, as done in 
the appendix, and to obtain the equation 

detM = (22) 

where the 3P x 3P hcrmitian matrix M is given by 

M„„,/j^ = {ojI - uj^)6a05^,,. + tJ^A„^,0^(O) (23) 

with a, /3 € {x, y, z}, e {1, • • • , P}, and lo^ is given 
by Eq. pSl) . We have defined 



2{u}b/cf Erfi(w&/c) 



KGRL 



27ri/2(^6/c)3 eMA) 
Sap A-'2 



KK'p 



(w/c)2-i^'2 



(24) 



where K' = K — q and the sum is over the reciprocal 
lattice of the underlying Bravais lattice. The index ^ or 
V labels the P atoms inside the primitive unit cell, and 
and are the positions of the corresponding atoms 
in that cell. By a slight generalization of the technique 



developed in [ij], it is possible to show that 



,,/3^(0) = /^afi,pu{il) e 



{u>b/c) 



(25) 



where Omin is the minimal distance between two atoms 
in the crystal. In particular, this implies that, for [i ^ 
the sum over K in p4l) has a finite limit when 6-^0. 

Here again, we note that frcc-Hcld solutions, that is 
eigenmodes of the free field with an electric field that 
vanishes in all atomic locations, may have been left out 
from (p^l) and should be investigated separately. The 
lower bound (121]) on the possible free-field eigenfrequen- 
cies Wfree Still applies [2^ . but, in the non-Bravais lattice 
case, this lower bound is in general not reached, since the 
vanishing of the electric field on each atomic location 
inside the primitive unit cell adds extra constraints with 
respect to the Bravais case. 



V. THE LONG WAVELENGTH LIMIT: 
GENERALIZED LORENTZ-LORENZ EQUATION 

Let us consider the long wavelength (LW) limit of the 
energy spectrum, for a general (even non-Bravais) lattice. 
We expand Eq. (l23l) in the limit ga <C 1 and wa/c <C 1, 
where a is the minimal distance between two atoms in 
the underlying Bravais lattice. The resulting equation is 



det] 



irLW 



0, 



(26) 



where the matrix 



(tLW 



itLW 



is given by 



at^,Pu - (t^o ^ ^ )SapS^i' 

, 2 SaP 



q»qp 



(lj/c)'^ — q^ 



(27) 



with 



KeRL* 



lim 

6-+0 



SapSfj.1, 



Vl 



127r3/2 63 

2 



Sgp K - KgKp _/f2j,2 
if2 ^ 



(28) 



The matrix Sapq — qaqp in Eq. ()27p is proportional to 
a projector: it has eigenvalues equals to zero or to g^, 
respectively for the modes which have a longitudinal or 
transverse polarization with respect to q. Equation (|26|) 
is a generalization of the result obtained by Hopficld in 
Eq.(30) of [l|, where the excitation spectrum in the long 
wavelength limit was calculated for a cubic crystal. Since 
the cubic crystal is a Bravais lattice, we can drop the 
indices fj, and v; using the symmetries of the cubic struc- 
ture one finds that Ja,p = 2(5^^/3. With this value one 
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obtains that Eq. ([26]) gives both the longitudinal mode 
frequency ignored in 



LW, cubic 

^11 ■ 



(29) 



and the transverse mode frequencies, solution as in [ij of 
the Lorentz-Lorenz (or Clausius-Mossotti) equation 



LW, cubic 
'± 



1 



/ LW. cubic \ 
nOL\ijJ ^ ■ j 



I \n 1 1 / LW.cubicx (^^) 

Icy 1 - 3 na(wj_ ' ) 



where n — 1/Vl is the atomic density of the cubic lattice 
and the real function a(w) is an atomic polarizabiUty 



nQ;((jj) 



(31) 



Then Eq. ()26p is a generalized Lorentz-Lorenz equation 
to any periodic structure. 

To illustrate these results, we have plotted in FigO 
the dispersion relation for light obtained in the long- 
wavelength approximation (|26p and from the exact equa- 
tion ([m, for a cubic lattice and for an orthorhombic 
lattice. As expected, the long wavelength approximation 
becomes inaccurate for 1 < qa. 

As a side remark, it is worth stressing that one can 
relate the tensor J to the dipolar electrostatic energy of 
the lattice. For a Bravais lattice L one obtains |23| 



1 



Jq/3 ^ 3^"^^ 



E 



(32) 



The sum in the right-hand side of this equation is not 
absolutely convergent, so that its precise meaning needs 
to be defined, as we did in I23l. 



VI. THE DIAMOND CRYSTAL: 
OMNIDIRECTIONAL PHOTONIC BAND-GAP 

As pointed out in the introduction, no atomic Bra- 
vais lattice leading to an omnidirectional gap for light 
is known. A natural idea to obtain such a gap is thus 
to consider lattices with several atoms per primitive unit 
cell. If one replaces the point-like atoms with extended 
objects like dielectric spheres, it is known that the dia- 
mond lattice leads to a photonic gap provided that the 
spheres fill a large enough fraction of the unit cell jlT| . 

We have therefore calculated the band structure for 
light in an diamond atomic lattice, solving numerically 
of section IIVI The diamond lattice is formed 
by the superposition of two copies of the same Bra- 
vais lattice: the fee lattice of lattice constant a, gen- 
erated by the three basis vectors ei = (0, a/2, a/2), 
©2 = (a/2,0, a/2), 63 = (a/2, a/2,0), and a second fee 
lattice obtained by translating the first lattice by the 
vector (a/4, a/4, a/4). For the reciprocal lattice of the 



fee lattice, we take the basis ei = (— 27r/a, 27r/a, 27r/a), 
§2 = (27r/a, — 27r/a, 27r/a), 63 — (27r/a, 27r/a, — 27r/a). In 
view of application to atomic gases, we restrict to the 
perturbative limit Up/ujQ — > 0, around the atomic reso- 
nance frequency ujq, calculating w — cjq to first order in 
ujp. In this regime, it is convenient to take as a unit of fre- 
quency the spontaneous emission rate F of a single atom, 
defined in (|20p . and one finds that lu — ujo is any of the 
eigenvalues of the matrix Vafi,pv{0) = linib^o IPa/j,/3i^(^)j 
with 



l + 2(fco6)' 



+ 



SttF 



KeKL 



27rV2(fco5)3 



- Erfi (fcofe)e 



(33) 



where K' = K — q, ko = luq j c is the resonant wavevec- 
tor, and where the limit for & ^ is computed in prac- 
tice from the relation PQ^,/3,y(0) = PQ/i,/3iy(6) e'^o'' for 
h ^ amin = aV3/4 (h = 0.05a in the present calcu- 
lations). Clearly, in this perturbative limit, [lo — LOf^jV 
only depends on the dimensionless quantities qa and k^a. 

We investigated first the arbitrarily chosen value k^a — 
2. In FigU^. we give w — £^s a function of q along the 
standard irreducible path in the first Brillouin zone of 
the fee lattice: a gap appears. To prove the existence 
of a truly omnidirectional gap, exploring the irreducible 
path is not sufficient. We thus evaluated the full density 
of states p(w), 



6 

E 



rf3 



9 



FBZ 



(27r)3 



6(UJ - Wq,„) 



(34) 



where the integral over q is taken in the first Brillouin 
zone, the sum runs over the band index n, tjq_„ is the 
light dispersion relation in band n. Numerically, we 
have explored the whole unit cell of the reciprocal lat- 
tice, parametrized as q = X^aQaea, Q G [—1/2,1/2)'^, 
with a cubic mesh for Q with 10^ points. The result in 
FigEl unambiguously shows the presence of an omnidi- 
rectional gap, piio) vanishing over an interval [winf, Wsup]- 



The values of 



(^inf 



and 



'^sup 



are reported as horizontal 



dashed lines in FigH^, showing that remarkably the par- 
tial gap on the irreducible path is very close to the full 
omnidirectional gap. 

In a second stage, we have varied the lattice spacing a. 



The variation of the gap A 



o^inf as a function 



of k^a is given in Fig|S^, revealing that there exists no 
gap for too high values of k^a^ 5.14 < k^a. This may be 
understood by looking again at the light dispersion rela- 
tion along the irreducible path of the first Brillouin zone, 
for values of k^a around 5.1, see FigdjD: Around point 
L, defined by = {-k / a,i: / a,TT / a), a frequency branch 
is seen to rapidly move towards low frequencies when a 
increases, leading to a closure of the gap. The modulus 
of the Bloch vector is qL = 7rv^/a ~ 5.44/a, and when 
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qL approaches the value fco from above for increasing a, positive frequency part 
the denominator — q\ appearing in ([55)1 for K = 



becomes very small and negative, leading to a mode fre- 
quency LJ close to the L point more and more below the 
atomic resonance frequency wq- In the opposite limit of 
a small lattice spacing a, that is of an increasing atomic 
density, one finds that the gap increases; from the overall 
factor l/(fcoa^) appearing in psp . one naively expects a 
gap scaling as Ija^ in this limit, which is indeed roughly 
the case. 

The previous discussion can be extended to the vari- 



ation of the gap boundaries Wjnf and 



^supj 



plotted as 



functions of fcoa in Figl6)3. The log-log scale reveals that 
LJini and Wsup approximately vary as power laws with a at 
low a, and the corresponding slopes ~ —3.25 and ~ —2.8 
(see the dotted lines) indeed indicate exponents close to 
—3. At large a, the rapid variation of ijJsup leading to the 
gap closure around Winf — = Wsup — ^^o — — T is also 
quite apparent. 

To be complete, let us mention that the previous dis- 
cussion about the existence of a spectral gap is not 
affected by the free-field solutions. From (PT|) we in- 
deed find that any free-field solution has a frequency 



> c infKGRL- K/2 = c|ei + 62 + e3|/2 = tt^/Sc/ 



a. 



We have seen that having a spectral gap requires fcg < 
— TrVs/a, which thus implies < uJbee |24l |. 



Experimental realizability of a diamond atomic 
lattice 



To experimentally test the presence of a gap for light 
in a diamond lattice with real atoms, the first step is 
to use an atomic species with a transition between a 
ground state with a spin Jg = and an excited state 
with a spin Jg = 1, so that the harmonic oscillator rep- 
resentation of the atomic dipole in the Hamiltonian H 
is reasonable. A natural candidate is strontium ^'^Sr, al- 
ready used in experiments on light coherent backscatter- 
ing |25| . but there are of course other possibilities, such 
as bosonic ytterbium ^^"'Yb where a Bose-Einstein con- 
densate is available [2^. The second step is to realize a 
diamond structure, by loading the atoms in an appropri- 
ate optical lattice with a filling factor one (Mott phase). 
Very recently, such a Mott phase with ^^^Yb was realized 
in a cubic lattice [l^l . 

Atomic lattices with two atoms per primitive unit cell 
have already been realized experimentally (28j . however 
to our knowledge not with a diamond structure. This 
missing element, the diamond lattice, can be realized 
building on the ideas of as was shown in [l^ in the 
different context of holographic lithography of dielectric 
materials [2^. We recall here briefly the idea: The op- 
tical lattice is produced by the coherent superposition of 
four laser plane waves, resulting in the total laser field of 



E 



(+) 

laser 



(r,t) = ^E,< 



,i(kj- 



(35) 



where the laser frequency tiiiascr is very far from the 
atomic resonance frequency coq so that the optical lattice 
is a purely conservative potential. The vectorial ampli- 
tudes of the laser plane waves have to satisfy the transver- 
sality conditions 



k, . E, = Wj. 



(36) 



The lightshift experienced by the atoms in this laser field 
is proportional to the laser intensity /lasoi (r) cx |E|j^^j.p. 
Taking for simplicity linearly polarized plane waves, the 
vectorial amplitudes Ej may be chosen real by a conve- 
nient choice of the origin of coordinates , so that 



IE 



(+) 12 

laser I 




J2 2E,.E,vcos[(k,-k,0T]. 

o<j<r<3 

(37) 

The laser wavevectors kj are chosen so that the kj — kj' 
generates the reciprocal lattice of the fee lattice. As in 
III one may choose ko = f (0,-2,-1), ki = f (2,0,1), 
k2 = f (0,2,-1), and kg = f (-2,0, 1). These four vec- 
tors have the same modulus, as it should be, which relates 
the value of the laser frequency to the diamond lattice 
constant. 



^laser 



\/57r 



(38) 



To have a significant spectral gap, see FiglSl one should 
satisfy the condition fega < 5, that is wiaser > 1 -41^01 
which corresponds to a blue detuned lattice, where the 
atoms are trapped in the minima of the laser intensity. 

The challenge is now to correctly choose the vectorial 
amplitudes E^ to ensure that the minima of /laser (r) form 
a diamond lattice. An elegant solution was given in [l|: 
by imposing the conditions 



3 • Ji/2 — rj3 • Jl/o = ■ii'2 • ■Cji 
E3 • El = E2 ■ Eq = 0, 

one realizes an intensity pattern 



-Ei-Eo>0, (39) 
(40) 



' laser 



(r) = /q + /i — cos(^ • r) + ^ cos(e^ ■ i 
7— 1 7— 1 

(41) 

where Iq > 0, /i > 0, and the basis vectors e.^ of the re- 
ciprocal lattice are given in the beginning of this section 
rVll The solutions of dMl) , dSS]) , (SOI) are not unique. If one 
fixes the value of Ji (that is the value of E3 • E2) to sup- 
press a global scaling invariance, we are left with 10 equa- 
tions for the 12 unknown components of the Ej's, leading 
to an actually explicitly calculable continuum of solu- 
tions parametrized by two real parameters. One of these 
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continuous parameters corresponds to the invariance of 
all the equations by the scaling transform E^- — > E^/A, 
j e {1,3}, and E^ ^ EjX, j e {0,2}, with the scaling 
factor A e R* . The particular solution given in O obeys 
the two constraints Ef = E2 = presumably added for 
experimental convenience. 

The key point is then that the minima of the inten- 
sity pattern (|1T|) are located in (3a/8, 3a/8, 3a/8) and 
in (5a/8, 5q/8, 5a/8) modulo any vector of the fee lat- 
tice [30|. These minima correspond to the same laser 
intensity, and lead to local harmonic microtraps that are 
isotropic. Since the relative vectorial position of these 
two minima is (a/4, a/4, a/4), the set of all laser inten- 
sity minima indeed form a diamond lattice. 



VII. CONCLUSION 

We have investigated the Fano-Hopfield model for 
propagation of light in a periodic ensemble of atoms with 
fixed positions for the most general (even non-Bravais) 
lattice geometry and without restricting to the long wave- 
length regime for light. We have shown that all diver- 
gences in the large momentum cut-off limit may be elim- 
inated in a systematic way by a renormalization of two 
parameters of the Hamiltonian, the atomic resonance fre- 
quency and the atomic dipole. 

First, we have used our theory to quantitatively con- 
firm the predictions of Q for the spectrum of light in a 
face-centered cubic lattice and the absence of an omni- 
directional spectral gap in that case. Second, we have 
shown that a diamond atomic lattice can lead to an om- 
nidirectional gap. Since a diamond atomic lattice may 
be realized in practice, as we discussed, this opens up 
the possibility of observing such a spectral gap for light 
in an ultracold atomic ensemble, which would bridge the 
state of the art gap between atomic gases and photonic 
band-gap solid state materials. 
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APPENDIX A: DERIVATION OF THE 
EXCITATION SPECTRUM IN TERMS OF BARE 
QUANTITIES 

1. Equations of motion 

Starting from H in and from the bosonic com- 
mutation relations it is possible to derive the system of 
coupled equations of motion in the Heisenberg picture 



for the matter excitation operators bi^a < i < N and 
a G {x, y, z}), and for the light excitation operators ak,e 
(k G X* and e _L k): 



(Al) 



~hwBbl^+M^A^^)^ (A2) 



ih — dke — hck ake 
at 

- Cu x(k; 6) E E (^^" + ^U) ' (A3) 

i—1 a 
N x,y,z 

- Cu x(k; E E + bl j e^^'^^ , (A4) 

where Ck = -dgSk and £k = i[hkc/{2eoV)]^/'^ . The 
standard Bogoliubov procedure to obtain the eigenmodes 
amounts to replacing the operators in the equations of 
motion with complex numbers oscillating in time at the 
eigenfrequency 



(A5) 
(A6) 
(A7) 
(A8) 



The set of equations (jAip - (jA4[) then becomes 



(A9) 
(AlO) 



N x,y,z 

- Ck x(k; 6) ^ £ ([/..„ + V,,a) el e-''^-^%(All) 

i—1 a. 

^'^ke = —hckvitie 

N x,y,z 

- Ck x(k; &) E E (f^^." + ^^.") ' 



where 



dBEi,„(R,) = 

(A13) 



keX'e-Lk 



Equations (|X9)) - (|M2)) can be rewritten, after simple 
manipulations, in the useful form 



, . l"^"" dB Ex,^(R,), (A14) 



OJB 



(C/,,„ -FV-,„), 



(A15) 
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AT x,y,z 



hick + Cij) 



Ea. (|A14|) is valid for 



^ 0, and Eq. fJIe)) is valid 



for ck — uj ^ 0. We note that a solution with a; = c/c is 
possible if one finds a free field solution (that is a solu- 
tion of Maxwell equations of frequency to in the absence 
of charge and current) such that the electric field Ej^ van- 
ishes in all atomic locations R,;, so that all Ui^a + Vi,a are 
zero. This type of free field solutions, already consid- 
ered in Q, are thus left out by Eq. (|A16[) and have to be 
investigated by a specific calculation. 



Inserting (jA20[) and (|A2ip in the expression of Ex^Q.(Ri) 
of Eq.jMal, using the relation X/£±k^«^^ — ^(^P — 
kakfj/k"^, and Ea. (|A14[) . and using the fact that x(k; h) is 
an even function of k, we obtain that the eigenfrequency 
ijj is given by the solution of the equation detM(6) = 
where the matrix elements Maj3{b) are given in Eq.®. 
3. Periodic system: non-Bravais lattice 



Finally, we briefly consider the case of a generic non- 
Bravais lattice of section IIVI in the elementary unit cell, 
there are P atoms of positions r^, v S {1, . . . , P}, and 
this base is repeated periodically according to an under- 
lying Bravais lattice. In the equations of subsection lA II 
one then has to replace the Bravais lattice positions 
by the atomic positions in the crystal: 



2. Periodic system: Bravais lattice 

We now consider TV atoms in a generic Bravais lat- 
tice, within the quantization volume V = {r | r = 
2:76^,0 < < M^}, 7 e {1,2,3}, e N*, com- 
patible with the lattice geometry. Thanks to the Bloch 
theorem one has 



i"j = r^. -f Ri 



(A22) 



where R^ belongs to the underlying Bravais lattice and 
is a position within the elementary unit cell. Then 
Bloch theorem gives 



(A18) 
(A19) 



where q G I? is the Bloch wavevector in the first Bril- 
louin zone in the reciprocal lattice (RL) and atom 
is placed at the origin of coordinates. By substituting 
(|A18[) and (|A19p in the expressions for Wke and Wke of 
Eas. (jA16[) and (|A17p . and using the fact that the sum 
'^j^^i g±«(k±q) Ri _ Q values of k± q that are not 

points of the reciprocal lattice, while ^^j^ g±«(k±q) Ri _ 
V/Vl for k ± q = K € RL, where V is the quantization 
volume occupied by the N atoms and Vl is the volume 
of the primitive unit cell in the direct lattice L, one has 
that Uke = and = except if k = q + K, K e RL, 
in which case 



(A23) 
(A24) 



where the Bloch vector q e P may be chosen in the first 
Brillouin zone of the reciprocal lattice RL. For k ~ q+K, 
where K € RL, 



Uk£ = 



Vl h{ck — uj) 



x^y.z P 



(A25) 

otherwise Uke is equal to zero. Similarly, for k = —q+K, 
where K e RL, 



or if k = — q + K, K G RL, in which case 



Ck X(k; b) 
^■"^ Vl h{ck + uj) 



(A20) 



(A21) 



_ Ck X(k; b) Y^'^rrjM ,y(-))^k-r 

""^^ " Vl Kck + c.) ^ ^ ^^"'^ + ^ ' 

(A26) 

otherwise Uke is equal to zero. Then proceeding as for 
the Bravais case one gets 



[1] J.J. Hopfield, Phys. Rev. 112, 1555 (1958). 

[2] V. Agranovich, Sov. Phys. JETP 37, 307 (1960). 

[3] U. Fano, Phys. Rev. 103, 1202 (1956). 

[4] P. de Vries, D.V. van Coevorden, A. Lagendijk, Rev. 

Mod. Phys. 70, 447 (1998). 
[5] J. Kastel, M. Fleischhauer, and G. Juzeliunas, Phys. Rev. 



A 76, 062509 (2007), and references therein. 
[6] I. Carusotto, M. Antezza, F. Bariani, S. De Liberato, and 

C. Ciuti Phys. Rev. A 77, 063621 (2008). 
[7] D.V. van Coevorden, R. Sprik, A. Tip, and A. Lagendijk, 

Phys. Rev. Lett. 77, 2412 (1996). 
[8] J. A. Klugkist, M. Mostovoy, and J. Knoester, Phys. Rev. 



9 



Lett. 96, 163903 (2006). 
[9] M. Greiner, O. Mandel, T. Esslinger, T.W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[10] K. I. Petsas, A. B. Coates, and G. Grynberg, Phys. Rev. 

A 50, 5173 (1994). 
[11] K.M. Ho, C.T. Chan, and CM. Soukoulis, Phys. Rev. 

Lett. 65, 3152 (1990). 
[12] J.D. Joannopoulos, S.G. Johnson, J.N. Winn, R.D. 

Meade, Photonic Crystals: Molding the Flow of Light, 

2nd ed. (Princeton Univ. Press, 2008). 
[13] E. Yablonovitch, T.J. Gmitter, K.M. Leung, Phys. Rev. 

Lett. 67, 2295 (1991). 
[14] M. Antezza and Y. Castin, arXiv:0903.0765 ^^1 (2009). 
[15] O. Toader, T.Y. Chan, and S. John, Phys. Rev. Lett. 92, 

043905 (2004). 

[16] The use of a Gaussian cut-off ensures a rapid convergence 
of the sums to appear over the reciprocal lattice. Com- 
bined with the relations (1191) and (|25|l . it allows to avoid 
Ewald summation tricks [ITI. [T^ . with which it shares 
however some mathematical features. 

[17] G.D. Mahan, J. Chem. Phys. 43, 1569 (1965). 

[18] M.R. Philpott, J.W. Lee, J. Chem. Phys. 58, 595 (1973). 

[19] Not all cut-off functions lead to an expression for (jSj 
which is possible to renormalize when 6 — > 0. Surpris- 
ingly, even by considering the simple case of a simple 
cubic lattice of lattice constant a, if one chooses for 
x(k; h) the function equal to one for k inside the cu- 
bic box [— f , f^]^ and zero elsewhere, a/h being inte- 
ger, one obtains that, when b ^ 0, the sum (|8]) di- 
verges, and, most important, the divergent terms de- 
pend on the Bloch wavevector q, hence providing a non- 
renormalizable model. To show this, let us consider the 
sum (|8} for a = x and f3 — y. By using symmetry prop- 
erties of the cubic geometry the sum has, for b ^ 0, a 
leading term 

KeRLn[-f,fl3 ^ ^ 

that diverges in a way proportional to q^qy for b ^ 0. 

[20] L.H. Ryder, in Quantum Field Theory (Cambridge Uni- 
versity Press, 1996). 

[21] For loq to be positive, b should be large enough. For an 
hydrogen- like atom we may use the estimates huiB ~ 
g^/(47reoao) and ds ~ \l\o-o, where q is the electron 
charge and ao the atomic radius, so that d%/{h£ociJB) ~ 
47raQ. The first factor in (|16p is positive if ao < 6, which is 
natural in a dipolar coupling model. Then the second fac- 
tor is close to unity if (ao/b) x (a;sao/c)^ < 1. This condi- 
tion is satisfied, since uosao/c ~ a, where a ~ 1/137 <^ 1 
is the fine structure constant. 



[22] The electric field of a free-field solution is a sum of plane 
waves which interfere destructively in all atomic posi- 
tions. This sum thus includes at least two plane waves. 
Let us select arbitrarily two of these plane waves; if one 
has a wavevector k, the other one has a wavevector nec- 
essarily of the form k + K, where K G RL*, because of 
the Bloch theorem. These two waves have the same fre- 
quency uJhcc and obey the free space dispersion relation, 
so that ojfrco/c = k = |k-|-K|. Squaring this last equality 
gives 2k • K -I- Jl^ = 0, so that k = kx - K/2 where kx 
is orthogonal to K. Then k > K/2, which leads to ([21}. 

[23] One introduces Sai3(k; b) — "^ZueL" exp(ik • IV){uap{r + 
R))r, with Uafsir) = {Sap - 3rar^/r^)/(47rr^) and {. . .)r 
denotes the average over r with the Gaussian proba- 
bility distribution oc exp(— r^/4b^). Using Poisson's for- 
mula and restricting to fc ^ iufKefli,* K, one obtains 
limt^o 5*0/3 (k; b) ~ - S^p + J^^g] /Vl. The Fourier 

k k J 

transform of Ua/s is indeed Hap — ^-i — ^Sap- We define 
the right-hand side of (|32p as limi,_o(Sct/3(k; b))k, where 
(. . .)k stands for the uniform average over the direction 
of k. We then get 

[24] A mathematical question is to know what is the ex- 
act minimal free-field frequency in the diamond atomic 
lattice. Within the two-mode ansatz of one finds 
a minimal ajfroo/c = 2Tv\/^/a. The general three-mode 
ansatz can give lower frequency solutions that however 
form a discrete set. The corresponding minimal value is 
i^froc/c — Stt/ (a\/2); it is obtained by superimposing with 
the right amplitudes the three plane waves of wavevec- 
tors k = (— 37r/2a, 37r/2a, 0), k — ei and k -I- §2, with a 
common polarization e = (1, 1,0). We have not explored 
the ansatz with four waves or more. 

[25] Y. Bidel, B. Klappauf, J.-C. Bernard, D. Delande, G. 
Labeyrie, C. Miniatura, D. Wilkowski, R. Kaiser, Phys. 
Rev. Lett. 88, 203902 (2002). 

[26] Y. Takasu, K. Maki, K. Komori, T. Takano, K. Honda, 
M. Kumakura, T. Yabuzaki, and Y. Takahashi, Phys. 
Rev. Lett. 91, 040404 (2003). 

[27] T. Fukuhara, S. Sugawa, M. Sugimoto, S. Taie, Y. Taka- 
hashi, Phys. Rev. A 79, 041604 (2009). 

[28] M. Anderlini, P.J. Lee, B.L. Brown, J. Sebby-Strabley, 
W.D. Phillips, J.V. Porto, Nature 448, 452 (2007). 

[29] A. Chelnokov, S. Rowson, J.-M. Lourtioz, V. Berger, J.- 
Y. Courtois, J. Opt. A: Pure Appl. Opt. 1, L3 (1999). 

[30] The locations of the intensity maxima are 
±(a/8, a/8, a/8) modulo any vector of the fee lat- 
tice, so that they also form a diamond lattice. 



\ ' / 

\ ' / 

""""^^^s^ 

1 ^ — 
1 


/ 1 \ ' V 

/ ' \ ' 
/ 1 \ 1 
/ ' \' 
/ ' \ 

/ ' / 
/ 1 /I 

1 /I 
\_ 1 / 1 

1 1 




1 

1 1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


\ 1 / 1 
\ ' /I 
\ ' / ' 
\ ' / 1 
\ 1 / 1 
\ 1 / 1 
\ 1 / 1 
\ ' / ' 
\ ' / ' 
\ 1 / 1 
\ 1 / 1 
\ 1 / 1 
\ 1 / 1 
\ 1 / 1 
\l/ 1 
V 1 






FIG. 1: Dispersion relation for light in a fee atomie lattiee, 
with fcoffl ~ 5, r/tjo ~ 0.0167, where a is the lattiee eonstant 
and r the spontaneous emission rate defined in H20p. In (a) 
the predietion found in figure (1) of [3, in (b) the results of 
calculations obtained using Eq. (|18|) of this paper, and (c) is a 
zoom of (b) around the atomie resonance frequency, (b) and 
(c) show the absence of an omnidirectional photonic gap, in 
contradiction with (a). 
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FIG. 2: (Color online) Dispersion relation for light in a fee 
atomic lattice, with koa ~ 3.14, T/ujo ~ 1.189, where a is the 
lattice constant, (a) is the prediction found in figure (2) of 
[circles correspond to the free modes, i.e. modes of the 
free field that vanish at all atomic positions] , (b) is the result 
of calculations using Ea. (|18|l of this paper [the red-dashed 
lines correspond to the free modes] . (a) and (b) appear to be 
identical. 
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FIG. 3: (Color online) Dispersion relation for light, for the ap- 
proximate long wavelength expression of Eq. (|26p (red), and 
for the complete expression Eq. (|18|) (black). The spectrum is 
calculated as a function of q/ko along direction (1, 1, 1) in re- 
ciprocal space. In (a) the simple cubic lattice with koa =1/5 
and V/ujo = 5 x 10"^, where a is the lattice constant. The 
long wavelength approximation fails when the horizontal (de- 
generate) branches of the exact spectrum start approaching 
each other, eventually closing the gap at larger values of q/ko. 
In (b) an orthorhombic lattice with basis vectors in the direct 
space along x, y, z with norms respectively proportional to 27r, 
TT, and 2^''^7r; here fcoa = 1/5 and F/cJo = 5 x 10~^, where 

1/3 

a = Vl and Vl is the volume of the primitive unit cell. Five 
non-degenerate branches are present, and the long wavelength 
approximation fails when the three horizontal branches of the 
exact spectrum start approaching each other, eventually clos- 
ing the gap at larger values of g/fco. 
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FIG. 4: (Color online) Dispersion relation for light in a dia- 
mond atomic lattice as a function of the Bloch vector along 
the standard irreducible path in the first Brillouin zone, in 
the perturbative regime tOp/ujQ ^ 1 (see text), (a) koa = 2, 
with ko = ujq/c and a is the diamond lattice constant (see 
text). Symbols: numerical values of the allowed frequencies 
uj. Dashed horizontal lines: lower and upper borders of the 
gap obtained from the full density of states (not restricting 
to the irreducible path), (b) For several values of koa close to 
the vanishing of the gap: koa — 5 (black symbols), koa — 5.1 
(red symbols), koa — 5.2 (green symbols). 
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FIG. 5: Density of states p{uj) for light in an diamond atomic 
lattice, in the perturbative regime ujp IujI < 1, see JMl)- The 
whole unit cell of the reciprocal lattice is explored to obtain 
this density of states, with a mesh of 10® points. The his- 
togram (with 250 bins) clearly reveals the existence of an om- 
nidirectional gap at frequencies below the atomic resonance 
frequency. We have taken fcoa = 2 as in Fig|3^. p(i>S) is in 
units of l/(rVL) where F is the single atom spontaneous emis- 
sion rate and Vl = a''/4 is the volume of the unit cell of the 
fee lattice. 
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FIG. 6: (Color online) For the light spectrum in a diamond 
atomic lattice, in the perturbative limit ujp ^ ujq, and as a 
function of koa, (a) value of the omnidirectional spectral gap 
A — LOsup — i^int, and (b) positions of the lower frequency 
tJinf and upper frequency uisup of the gap with respect to the 
atomic resonance frequency, a is the lattice constant (see 
text), fco — ljq/c and F is the spontaneous emission rate de- 
fined in (|20p . The inset in (a) is a magnification. 



